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ABSTRACT 


Context. The solar brightness varies on timescales from minutes to decades. A clear identification of the physical processes behind 
such variations is needed for developing and improving physics-based models of solar brightness variability and reconstructing solar 
brightness in the past. This is, in turn, important for better understanding the solar-terrestrial and solar-stellar connections. 

Aims. We estimate the relative contributions of the continuum, molecular, and atomic lines to the solar brightness variations on 
different timescales. 

Methods. Our approach is based on the assumption that variability of the solar brightness on timescales greater than a day is driven by 
the evolution of the solar surface magnetic field. We calculated the solar brightness variations employing the solar disc area coverage 
of magnetic features deduced from the MDI/SOHO observations. The brightness contrasts of magnetic features relative to the quiet 
Sun were calculated with a non-LTE radiative transfer code as functions of disc position and wavelength. By consecutive elimination 
of molecular and atomic lines from the radiative transfer calculations, we assessed the role of these lines in producing solar brightness 
variability. 

Results. We show that the variations in Eraunhofer lines define the amplitude of the solar brightness variability on timescales greater 
than a day and even the phase of the total solar irradiance variability over the 11-year cycle. We also demonstrate that molecular lines 
make substantial contribution to solar brightness variability on the 11-year activity cycle and centennial timescales. In particular, our 
model indicates that roughly a quarter of the total solar irradiance variability over the 11-year cycle originates in molecular lines. 
The maximum of the absolute spectral brightness variability on timescales greater than a day is associated with the CN violet system 
between 380 and 390 nm. 

Key words. Sun: activity — Sun: solar-terrestrial relations — Sun: magnetic fields — Sun: faculae, plages — Sun: sunspots — line: 
formation 


1. Introduction 

Regular spaceborne measurements of the total solar irradiance 
(TSI, which is the spectrally integrated solar radiative flux at one 
au from the Sun) started in 1 978 with the launch of the NIMBUS 
7 mission (iHovt et al.lll992h. They re vealed that the TSI varies 
on multiple timescales (IWillson et aI.l[T98lUFr6hlichll2005h . The 
most striking features of the TSI records are the ~ 0.1% mod¬ 
ulation over the course of the 11-year solar activity cycle and 
more irregular variations on the timescale of solar rotation with 
amplitudes of up to 0.3%. 

Measurements of the spectral solar irradiance (SSI, which 
is the solar radiative flux per unit wavelength at one au from 
the Sun) showed the variability on the same two timescales (see 
iFlovd et all IZOOl [Harder et al.l l2009t iDeland & Cebulal l2(m 
and references therein). They also indicated that the amplitude 
of the relative SSI variations strongly depends on the wavelength 
and generally increases towards short wavelengths. For exam¬ 
ple, the variability of the solar Ly-a line (121.5 n m) irradiance 
over t he course of the 11-year cycle reaches 100% ([Woods et al.l 
[ 2 OOOI) : i.e., it is three orders of magnitude higher than the TSI 
variability. 
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It is widely accepted now that the SSI variability on 
timescales of a day and longer is caused by solar surface 
magnetic fie l ds (iDomingo et al. [2009]; Krivoya_et alj j^l lal : 
lErinolli et al.l [20131; Isolanki et al.l 12013 ; Yeo et al.l [2014al) . H- 
though several other physical mechanisms for SS I variability 
have been proposed (see e.g. [Kuhn & SteinlllQ^ iKuhn et ^ 

mn). 


Magnetic fields emerge on the solar surface in the form 
of magnetic concentrations desc ribed well by flux tubes (see 
[Solankiiri993l : [Solanki et al.[[2006l for reviews). Large flux tubes 
usually form dark sunspots, while small flux tubes appear as 
bright features ensembles of which compose faculae and the 
network, while individual small flux tubes are seen as bright 
points when observed at high resolution. The effects from the 
dark and bright features on solar irradiance do not compensate 
each other and the imbalance depends on the solar disc coverage 
by magnetic features of different sizes and on the wavelength 
observed. The coverage changes with time, leading to the vari¬ 
ability of solar irradiance on timescales from days to decades. 
In recent years a number of models of solar irradiance variabil¬ 
ity have bee n created that are to some extent based on such a 


concent Ce.g. Fligge et al.[ 200(1 Krivova et alJ 20031; Lean et al. 

2005 

; Shanrro et al. 2011b; Fontenla et al. 1201 It [Bolduc et al. 

2014 

;lBall et alJl2014t Yeo et al.ll2014bh. 
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The interest in solar irradiance variability is not limited to 
the solar community. The terrestrial climate responds to the 
decadal variations in solar irradiance (see e . g. reviews bvjHaigri 
l2007HGrav eTZIbOloHErmolli et alJElTlLSolanki erai2013l) 
and there is also evidence for a longer term influence of so¬ 
lar activity on climate (e.g. iBond et ^1200Ih . Studies of solar 
irradiance are also of high importance for stellar astronomers, 
who have been comparin g it with the var i ability of other 
lower main sequence stars (iRadick et alJll99§: iLockwood et al.l 
l2007t[Hall et ai.|[2009l) . The interest in s olar-stellar comparisons 
has recently bee n rekindled (see e.g. iMcOuillan et al.l 1201 2t 
iBasri et ^ 12013h by the unprecedented precision of broad- 
band stellar photome try achieve d with the launch of the Kepler 
( Borucki et al.ll201^ and Corot (iBorde et al.ll2003t iBaglin et al.l 
120061) space mi ssions and the an ticipation of the upcoming 
PLATO mission (iRauer et al.ll2014l) . 

Changes in the solar spectrum can be decomposed into the 
variations in the continuum and in spectral lines. The spectral 
lines are present all over the solar spectrum and they are respon¬ 
sible for its extremely rich and sophisticated structure. In the 
EUV (10-100 nm) and far UV (100-200 nm) below 170 nm the 
solar spectrum contains strong emission lines, while longward 
of 170 nm its profile is determined by many millions of atomic 
and molecular absorption lines, also called Fraunhofer lines. It 
is hardly possible to find a spectral interval without Fraunhofer 
lines in the visible and near-infrared parts of the solar spectrum. 
In the UV the immense number of spectrally unresolved lines 
form the UV line haze which completely blocks the cont inuum 
photons (e.g. lCollet et al.ll2005HShort & HauschilddiTOOSh . 

The variations in the Fraunhofer lines in the full-disc 
solar brightness spectrum over the course of the activ- 
ity cycle were first revealed in the pioneering studies by 
iLivingston & Holwegerl(ll982h who analysed the Kitt Peak spec¬ 
tral records for seven Fraunhofer lines and found a decrease 
of equivalent width ranging from no variation (Si I 10827.1 


A line) to 2.3% (C I 5380 3 A line) from 1976 to 1980. 
ILivingston & Holwegerl (Il982h suggested that this decrease may 
be attribute d to a change in the ph otosph eric temperature gradi¬ 
ent (se e also lHolweger et aDl983h . Later lMitchell & LivingstonI 
(Il99ll) employed Fourier transform spectrometer (FTS) observa¬ 
tions at the McMath Telescope to show that the absorption lines 
between 500 to 560 nm exhibited a 1.4% decrease in ampli¬ 
tude and 0.8% decrease in equivalent width, at maximum com¬ 
pared with minimum of cycle 21. Interestingly. [Livingston et aP 
(l2007h could not confirm this finding in the next solar cycle for 
the unblended strong Fe I and the Na ID lines. This is, probably, 
a good indication of how challenging such obs ervations are (see 
also d iscussion in iPenza et al.ll 206 ^ . Instead. FLivingston et al.l 
(l2007h reported cyclic variation in the Mn I 539.4 nm line and 
the CN 388.3 nm bandhead. 

Efforts have also been made to estimate the effect of the 
Fra unhofer lines on the SSI variability with models. For exam¬ 
ple, [Unruh|et^ ([T99^ used spectra computed using the proba¬ 
bility distribution functions and an estimate of spot and facular 
disc area coverages during the solar maximum to show that spec¬ 
tral lines are the major contributor to the TSI variations in the 11- 
year activity cycle timescale rathe r than the continuum . In this 
study we follow up on the idea of lUnruh et alJ d 19991) and em- 
plo v the SATIRE (Spectral And Total Irradi ance Reconstruction, 
see iFligge et alJUOOOl: iKrivova et al.ll200^ model to assess the 
relative role of continuum and spectral lines contributions to the 
SSI variability. The availability of state-of-the-art measurements 
of solar disc area coverage of magnetic features and radiative 
transfer calculations of their spectra, that now account for the 


effects of non Local Thermodynamical Equilibrium (non-LTE), 
allows us to define the contribution of spectral lines to the irra¬ 
diance variability over the entire solar spectrum (including the 
UV) and on different timescales of variability. 

In Sect.l^we describe our model. In Sect.[2we show how the 
Fraunhofer lines affect the brightness contrasts between the dif¬ 
ferent magnetic features and quiet regions on the solar surface. 
In Sect. |4] we define the spectral profiles of the irradiance vari¬ 
ability on the solar rotational, 11-year activity cycle, and centen¬ 
nial timescales. In Sect.|5]we show how these three profiles look 
if we first exclude molecular lines from our calculations and then 
all Fraunhofer lines. The main results are summarised in Sect.j^ 


2. The model 


We base our calculations on the SATIRE model. Its branch 
SATIRE-S (with “ S” standing for the satellite era, see 
IKrivova et ^ 2011 hi) has been updated and refined over the re¬ 
cent years ' ilR^et~^.ll2n~n 1201 4 et al.ll2ni4hl) and it cur¬ 
rently replicates over 92% of the observed TSI variability over 
the entire period of spaceborne observations. 

In SATIRE magnetic features observed on the solar disc are 
divided into three classes: sunspot umbrae, sunspot penumbrae, 
and the combined faculae and the network. The part of the so¬ 
lar disc not covered by these magnetic features is attributed to 
the quiet Sun. The brightness of each component (i.e. one of the 
three magnetic feature classes or the quiet Sun) is assumed to 
be time-invariant, but depends on the wavelength and position 
on the visible solar disc. The solar irradiance is calculated by 
weighting the spectra of the individual components with corre¬ 
sponding disc area coverages, i.e.: 


S{t,A) = 


?// 

solar disc 


Ik (A, r)ak(t,r)d£l. 


( 1 ) 


where the summation is done over the three SATIRE classes of 
magnetic features and the quiet Sun. For each of the compo¬ 
nents the integration is performed over the visible solar disc. 
Here Ik (T, r) is the emergent intensity from the component k 
at the wavelength A along the direction r. The functions ak(t, r) 
are the fractional coverages of the solar disc by the component 
k along the direction r, so that ak{t, r) dQ. gives the elementary 
solid angle covered by the component k. 

The two main building blocks of our model are the spec¬ 
tra Ik (T, r) and the fractional coverages ak{t, r) of the quiet Sun 
and mag netic featu r es. W e utilise the fractional coverages de¬ 
duced bv iBall et alJ d2012l) from the full-disc continuum images 
and magnetograms obtained by the Michelson Doppler Imager 
onboard the Solar a n d Heliospheric O bservatory (SOHO/MDI; 
IScherrer et al.lll995h . iBall et all (1201 2l) successfully used these 
fractional coverages to reproduce 96% of TSI variat ions dur¬ 
ing cy cle 23 (rep resented by the PMOD composite by iFrohlichl 
I2OO6I) . Following iBall et ^(1201 lLl2012l) we replace the integra¬ 
tion in Eq. ([1]) with the summation over magnetogram pixels. 
Then the ak{t,r) is substituted by akij{t) which represents the 
coverage of the pixel with abscissa i and ordinate j by the com¬ 
ponent k.The akij(t) values lie between 0 and 1 for the quiet Sun 
and faculae, and either 1 or 0 for umbra and penumbra. 

The brightness spectra of magnetic features and the quiet 
Sun 4 (A, r) are calculated with Non local thermodynamic 
Equilibrium Spectral SYnthesis Code (NESSY, Tagirov 2015, 
private comm.), which is a further developme nt of the COde 
for Solar Irradiance (COSI. IShapiro et^l2010() . This is in con¬ 
trast to previous SATIRE publications, which have been based 
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on s pectra calculated w i th the ATLAS9 code by iKuruczl (Il992h 
and ICastelli & Kurucj (Il994ll . NESSY simultaneously solves 
the statistical equilibrium equations for the elements from hy¬ 
drogen to zinc, taking the coupling between level populations of 
different elements via the electron concentration and radiation 
field into account. The employment of NESSY allows us to re¬ 
calculate the 4 (A, r) spectra over the entire frequency domain 
(in particular, including the UV, see Sect. IH with and without 
molecular/Fraunhofer lines. 

We adopt the ID temperature and density structures of the 
quiet Sun, faculae, s unspot umbra, and sun spot penumbra from 
Models C, P, S by iFontenla et ^ (Il999h . and Model R by 
iFontenla et al.l (l2006l) . respectively. The spectra emergent from 
magnetic features are calculated neglecting the Zeeman split¬ 
ting. The populations of atomic and molecular levels, as well as 
the electron concentration are self-consistently calculated with 
NESSY. The atomic linelist is compiled from the “long” linelist 
provided by Kurucz (2006, pr ivate comm.) and the V ienna 
Atomic Line Database (VALD. [Kupka et^ll999L 2000 1. The 
molecular linelist is complied from the list by kuriiczldlOOl 
and the Solar Radiation Phy sical Modeling database (SRPM, 


_ hy si 

see e.g. iFontenla et al.l 1201 ih . The strongest vibrational bands 
of the CN violet system and CH G-ba n d are calculated usin g 
the molecular consta nts bv lKrupp|(ll974h: Knowles et al.l(ll988h : 
[Wallace et al.l d 19991) (see also Shapiro et al. 2010 , 201 lai and 
references therein for more details). 

The full disk brightness spectrum of the quiet Sun calculated 
with such a setup is in good agreement with various spaceborne 


measurements of the solar spectrum around activity minimum 
conditions (see IShapiro et alj|2010t iThuillier et al.ll20Tll |2014 
l2015[ for a detailed intercomparison). The centre-to-limb varia¬ 
tions in the quiet Sun brightnes s agree well with measurements 
in the UV dShapiro et alil2013l) and visible (Tagirov 2015, pri¬ 
vate comm.) spectral domains. For the experiments presented 
below we use NESSY spectra calculated with the full molecular 
and atomic linelist, spectra produced excluding molecular lines, 
and finally spectra calculated excluding all Fraunhofer lines, i.e. 
when the total opacity is given only by the free-free and bound- 
free processes. 


3. Effect of the Fraunhofer lines on contrasts of 
solar magnetic features 

In the upper left panel of Fig. [T] we plot the spectral fluxes at 1 
AU (hereafter referred as brightness spectra or just brightness) 
emergent from the quiet Sun, faculae, spot umbrae and penum- 
brae calculated as discussed in Sect.|3 For illustration purposes 
the brightness of each type of magnetic feature is plotted assum¬ 
ing that it fully covers the visible solar disc. 

With the exception of the umbral brightness, the spectral 
profiles of the brightness of different magnetic components are 
strongly affected by numerous absorption features. The umbral 
brightness is significantly lower than brightnesses of all other 
magnetic components considered in this study and the spectral 
lines have only marginal contribution to the umbral brightness 
contrast with respect to the quiet Sun. Hence the effect of spec¬ 
tral lines on the umbral spectrum will not be discussed in detail, 
though it will be taken into account when computing irradiance 
variability. 

In the lower left panel of Fig. [T] we present facular and 
penumbral brightness contrasts with respect to the quiet Sun. 
They are calculated by subtracting the brightness of the quiet 
Sun from the facular and umbral brightness. The overall profiles 


of the contrasts are to some extent reminiscent of the solar spec¬ 
trum but skewed towards the UV because at short wavelengths 
(i.e. when the Planck function can be approximated by the Wien 
approximation) the derivative of the Planck function with respect 
to the temperature (which defines the contrast) is proportional to 
the Planck function divided by the wavelength. 

To illustrate the contribution of spectral lines to the facu¬ 
lar and penumbral brightness contrasts we also show brightness 
spectra and contrasts calculated by putting the opacity of molec¬ 
ular lines to zero (middle panels of Fig. [TJ as well as putting 
opacity in all Fraunhofer lines to zero (i.e. producing purely con¬ 
tinuum spectra and contrasts, see right panels of Fig. [1]). We dis¬ 
cuss these plots in Sects. IXTI and B.21 respectively. 


3.1. Molecular lines 

The middle panels of Fig. [T] illustrate that the elimination of 
molecular lines from the radiative transfer calculations has a 
strong effect on the brightness spectra and the contrasts in the 
300 - 450 nm spectral domain. The strongest peak in the facular 
contrast is linked to the CN violet system (B^E - transi¬ 
tions) between 380 nm and 390 nm, while the second strongest 
peak results from the CH G-band (A^A-X^Tl transitions) around 
430 nm. These molecular systems, being also prominent features 
of the solar spectrum, almost double the facular contrast at the 
corresponding wavelengths. The brightness contrasts and spec¬ 
tra between 300 and 350 nm are also noticeably affected by the 
CN, NH, and OH bands. 

The amplification of facular contrast by molecular lines is 
due to the strong sensitivity of molecular concentrations to tem¬ 
perature changes. Molecular concentrations are lower in facu¬ 
lae which are generally warmer than the quiet Sun surround¬ 
ings. Consequently, the molecular lines in the facular spectrum 
are weaker than in the spectrum of the quiet Sun, which leads 
to the enh ancement of the contrast. Such an interpretation i s 
in line with lSteiner et al.l(l200ll) : ISanchez Almeida et al.l(l200ll) . 
Similarly to this study, they employed ID models of faculae 
and quiet Sun to explain the appearance of the G-band bright 
points, observed on G-band filtergrams, by the weakening of CH 
lines in hot flux tubes . Similar conclusion was later drawn by 
ISchiissler et al.l (l2003l) who employed more r ealistic 3D MHD 
simulations (see also iRiethmilller et aHl20l4 who studied the 
CN bright points). 

The CN violet system and CH G-band induce a compara¬ 
ble amplification of the facular brightness contrast even though 
the CH G-band is more pronounced in the solar spectrum. This 
is partly caused by the relatively high dissociation energy of 
the CN molecule (7.72 eV), which is larger than that of CH 
(4.25 eV). As a consequence the CN concentration is more sen¬ 
sitive to temperature changes than the CH concentration (see 
e.g. the discussion of t he chemical equilibrium calculations in 
iBerdvugina et al.ll2003l) . 

The CN violet system and CH G-band are less pronounced 
in the penumbral spectrum than in the spectrum of the quiet 
Sun, which is probably associ ated with a sm a ll grad ient of pho- 
tospheric temperature in the IFontenla et al.l (l2006l) penumbral 
model (see their Fig. 6). As a result the penumbral brightness 
contrast drops in the CH G-band and is only marginally affected 
by the CN violet system. 
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FULL LINELIST NO MOLECULAR LINES 


NO LINES 






Wavelength [nm] Wavelength [nm] Wavelength [nm] 


Fig. 1. Upper panels: the brightness spectra of the quiet Sun (black curves), faculae (red curves), spot penumbra (blue curves) and 
umbra (magenta curve, only plotted in the upper left panel). Lower panels: the facular (red curves) and penumbral (blue curves) 
brightness contrasts with respect to the quiet Sun. The dotted curves are calculated taking the spectral lines (full linelist) into account. 
The solid curves are calculated by putting the opacity in the molecular (middle panels) and in all spectral lines (right panels) to zero. 
To facilitate the comparison between different spectra and contrasts the dotted curves are plotted in all three panels. The two vertical 
dotted lines in the left panel constrain the spectral interval shown in the middle panels. 


3.2. Fraunhofer lines 

The right panels of Fig. [T] illustrate that the elimination of all 
Fraunhofer lines from the radiative transfer calculations has a 
profound effect on both the brightness spectra and the contrasts 
below about 600 nm. The spectral lines are a dominant source of 
opacity in the UV and they significantly contribute to the total 
opacity in the visible spectral domain. By shifting the forma¬ 
tion height of the radiation to higher and colder photospheric 
layers, the Fraunhofer lines significantly reduce the brightness 
of the quiet Sun, faculae, and penumbrae (see upper right panel 
of Fig.[T]l- The salient features of the continuum spectra are the 
Balmer, Mg I, and All ionisation edges at 364.5 nm, 251.2 nm, 
and 207.8 nm, respectively. The fe ature around 290 nm is at¬ 
tributed by iHaberreiter et al.l (l2003l) to a resonance in the Mg I 
3 photoionisation cross section. While the Mg I and Al I edges 
can still be recognised in the “real” (i.e. calculated with the full 
linelist, see dotted lines in Fig. [TJ solar spectrum, the Balmer 
edge is completely hidden in the multitude of spectral lines. 

The difference between the quiet sun and faculae continua is 
small only, the former slightly brighter in the middle and near 
UV (200-300 nm and 300^00 nm, respectively), the latter 
slightly brighter shortward of 200 nm and longward of 400 nm. 
The penumbra continuum is well below the other two, as can be 
expected from its effective temperature that is about 250K lower 
than that of the quiet sun. 

The Fraunhofer lines substantially increase the facular con¬ 
trast and decrease the absolute value of the penumbral con¬ 
trast. This is linked to the temperature structure of the ID at¬ 
mospheric models (see Sect. |2]i: the temperature difference be¬ 


tween the quiet Sun and faculae decreases, while the temperature 
difference between the quiet Sun and penumbr a increases, to¬ 
wards deeper photospheric layers (see Fig. 6 from lFontenla et aP 
l2006h . where the continuum intensity is formed. The effect of 
Fraunhofer lines on the brightness spectra and contrasts gets 
smaller in the infrared where the relative contribution of lines 
to the total opacity diminishes. 

According to our calculations the facular contrast in the con¬ 
tinuum becomes negative in the near and middle UV. The con¬ 
tinuum irradiance at these wav elengths com es from the deep¬ 
est photospheric layers (see also lAvrelll989l) . Indeed, at shorter 
wavelengths the photoionisation opacity from metals increases, 
while at longer wavelengths the opacity fr om negat i ve hy - 
drogen ions get higher (see Fig. 4.2 from iMihalasI Il978h . 
In deep photospheric layers faculae are cooler than the quiet 
Sun according to the ID atmospheric structures employed here, 
which results in the inversion of the contrast. A similar in- 
version in the continuum co ntrast was recently obtained by 
ICriscuoli & UitenbroekI d2014l) who based their calculati ons on 
3D MHD simula tions calculated by Fab bian et al.l (1201 2h using 
STAGGER code (iGalsgaard & Nordlundll996lf 


The negative full-disc facular continuum brightness contrasts 
in the UV obtained from our calculations are qualitatively con¬ 
sistent with the observations of negative continuum brightness 
contrasts of fac ulae about disc c entre in the visible part of the 
solar spectrum (lYeo et al.ll2013L and references therein) where 
the radiation also emanates from the dee p photosphe re. We note 
that the contribution functions (see e.g. lGravlll992L p.l51) for 
the full-disk UV and disc centre visible continuum radiation are 
comparable because the formation height of the continuum ra- 
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TIMESCALE OE SOLAR ROTATION 



ACTIVITY CYCLE TIMESCALE 



CENTENNIAL TIMESCALE 



Wavelength [nm] 


Fig. 2. The spectral profiles of the irradiance variability calcu¬ 
lated on the solar rotational (upper panel), activity cycle (mid¬ 
dle panel), and centennial (lower panel) timescales. Plotted are 
the values calculated employing “SATIRE-ATLAS9” SSI and 
TSI time series (blue curves) and values calculated employing 
the “SATIRE-NESSY” SSI and TSI time series (red curves). 
All spectral profiles are smoothed by applying a 10 nm running 
mean. 


diation decreases with cosine of the heliocentric angle but in¬ 
creases from the UV to visible. At the same time we are aware 
that while ID models are designed to repli cate many features 
of brightness spectr a and contrasts (see e.g. lUnruh et anil999t 
iFontenla et al.l20lTh they miss th e real physics of centre-to-limb 
variations in the facular contrast (ISteineij|2005l) . 

4. Spectral irradiance variability on different 
timescales 

The most accurate and uninterrupted records of solar irradi¬ 
ance variability available so far a re the TSI records (see e.g. 
the most recent publications by IKopd et ^ 12012l iFrohlichl 


1201f^hmutz et^l2013l:lKopt)ll20T4[ and references therein). 
Nevertheless more information than just TSI is needed for both 
climate simulations and the solar-stellar comparison. Climate 
simulations are very sensitive to the variability of the UV irra¬ 
diance, which is the main driver in the top-down mechani sm for 
the influence of the solar irradi ance variability on climate ( jHai iE 
Il994l^dera & Kurodal2003l . Likewise, most stellar photomet¬ 
ric records used for solar-stellar comparison are obtained with 
broad-band biters sensitive to the visible part of the spe ctrum 
(iLockwood et al.ll^OOTl : iHall et al.ll200^lBasri et al.ll2013l) . It is, 
therefore, important to know the spectral profiles of the solar ir¬ 
radiance variability on different timescales. In this Section we 
explain how these spectral profiles are linked to the brightness 
contrast discussed in Sect. |3] Then in Sect. |5]we show how these 
profiles are affected by the Fraunhofer lines. 

We define the spectral profile of the irradiance variability on 
the timescale of solar rotation as: 

trSSI ^ rms(SSI(T,0 -<SSI(T,0)8 i) 
crTSr rms(TSI(f) -<TSI(0)si) ’ 

where “rms” stands for the root mean square. In our study both 
TSI and SSI time series are calculated with daily cadence so that 
the averaging in Eq.|2]represents the running mean over 81 days. 

In the upper panel of Fig. we plot crSSI/crTSI values as 
functions of wav elength calculated with S ATIRE-S TSI and SSI 
time series from iBall et af] (1201 2l l2014l) . hereafter “SATIRE- 
ATLAS9” time series, since they are based on the spectra of in¬ 
dividual magnetic features and the quiet Sun calculated with the 
ATLAS9 code. These time series cover the period from February 
19, 1999 till October 2, 2010 and are based on the fractional 
coverages of magnetic features deduced from the SOHO/MDI 
data. We also plot crSSI/crTSI values calculated with the model 
described in Sect. |2] i.e. employing the same SOHO/MDI frac- 
tional coverages of magnetic features as iBall et~^ (I2012ll20l4) 
but convolving them with the NESSY spectra instead of the 
ATLAS9 spectra (see Sect. |2ll. These “SATIRE-NESSY” time 
series are calculated for the same period as the “SATIRE- 
ATLAS9” time series. 

In the middle panel of Fig.|2]we plot “SATIRE-NESSY” and 
“SATIRE-ATLAS9” spectral profiles of the irradiance variabil¬ 
ity on the 11-year activity cycle timescale. They are defined as 

ASSI ^ (SSI(T,0)2oo2-(SSI(T,0)2oo 8 
ATSr <TSI(f )>2002 - <TSI(f)> 2 oos ’ 

where annual averaging is performed over the calendar years of 
high (year 2002) and low (year 2008) solar activity. 

The irradiance variability on a centennial timescale is of¬ 
ten described as a secular change between solar min ima con¬ 
dition s (see e.g. the discussion of TSI variability in IFrohlichl 
l2009l) . The current solar irradiance records are too short and un¬ 
certain to unambiguously reveal and quantify secular changes. 
Consequently, the magnitude and even specific physical mech¬ 
anisms responsib le for the cen t ennial SSI variability are heav¬ 
ily debated (see ISolanki et al.l 1201 3l and references therein). 
According to SATIRE most of the SSI changes between activity 
minima are caused by the varyin g contribution from the network 
component (iKrivova et al.ll20lE) . The contribution of this com¬ 
ponent does not drop to zero even at activity cycle minima and is 
responsible for the secular trend in solar irradiance, e.g. for the 
change in irradiance between the 2008 activity minimum and the 
Maunder minimum (however, see also lSchriiver et all201 ll) . 
























































6 


A.I. Shapiro et al.: The role of the Fraunhofer lines in solar brightness variability 


Along these lines we define the spectral profile of the irradi- 
ance variability on the centennial timescale as 

ASSl _ {SSI(/l, 0)2008 “ SSIquiet(d) 

^ ^ “ (TSI(0>2008 - TSIquiet ’ 

where SSIquiet(d) and TSIquiet are spectral and total solar irra- 
diance, respectively, as they would be measured if the visible 
part of the solar disc did not contain any active features (i.e. was 
completely covered by the quiet Sun). We emphasise that our 
dehnition of the centennial variability refers to the irradiance 
changes between solar activity minima rather than averaged over 
the activity cycle as often done in the literature. The spectral pro- 
hles on the centennial timescale calculated with the “SATIRE- 
NESSY” and “SAT1RE-ATLAS9” time series are plotted in the 
lowest panel of Eig.|3 

In addition to strong concentrated magnetic helds, consid¬ 
ered in SATIRE, there is also wea k turbulent magnetic held on 
the so lar surface (see reviews by Ide Wiin et al.l 1200^ ISten 13 
I 2 OI 3 I) . Its effect on the rotational and 11-year irradia nce vari¬ 
ability is believed to be small (cf. ISolanki et al.ll201^ but it is 
presently unclear whether this held can contribute to the irradi- 
an ce variability on t he centennial timescale (see e.g. discussion 
in ijudge et al.ll20O and references therein). The ambiguity as¬ 
sociated with this effect is an additional source of uncertainty in 
our estimate of the contributions of molecular and atomic lines 
to irradiance variability on the centennial timescale. 

We note that the ATLAS9 code is based on the LTE as¬ 
sumption, which is not applicable in the UV. As a consequence 
iBall et af] (l2014t) corrected the SSI values between 115 and 270 
nm based on the empirical method bv iKrivova et al.l (l2006h . In 
contrast to ATLAS9, NESSY takes non-LTE effects into account 
and thus no additional correction of the “SATIRE-NESSY” time 
series is needed. Despite this difference in methods, the spec¬ 
tral prohles calculated with both time series are very close to 
each other over the entire spectral range shown in Eig. |2] also 
includin g the UV. We constru e this as an independent support to 
both the IKrivova et alJ (l2006h method and the reliability of the 
“SATIRE-NESSY” time series. 

The “SATIRE-ATLAS9” SSI and TSI time series have been 
demonstrate d to be consisten t with observation s from multiple 
sources (see iBall et al.l 120141 lYeo et alJl2014R and references 
therein). Since the differences between “SATIRE-NESSY” and 
“SATIRE-ATLAS9” spectral prohles are not essential for the 
purposes of this study we refrain from comparing the “SATIRE- 
NESSY” time series to observational data here, which is beyond 
the scope of this paper. 

Eigure |2 indicates that the spectral prohle of the irradiance 
variability depends signihcantly on the considered timescale. 
The SSI variability on the 11-year activity cycle timescale (mid¬ 
dle panel of Eig. |2]i is brought about by the competition between 
bright faculae and network on the one hand and dark sunspots 
on the other hand. The 11-year TSI variability as well as the 11- 
year variability in the UV and visible are faculae-d ominated and 
are in phase with the solar cycle (however, see also lHarder et alJ 
l2009h . so that the ASSI/ATSI values plotted in the middle panel 
of Eig. 13 are positive with the exception of a minor feature 
around 470 nm. The facular contribution decreases with wave¬ 
lengths faster than the spot contribution (in particular, the umbral 
component, see left upper panel of Eig. [T]| so that starting from 
600 nm both components almost cancel each other (while start¬ 
ing from 1-2 yum spots overweight the faculae a nd SSI varies 
out-o f-phase with solar cycle, see e.g. Eig. 7 from lErmolli et al] 
l2013h . As a result the dominant part of the irradiance variability 


over the 11-year activity cycle originates in the 250 - 450 nm 
spectral domain, which is especially affected by the molecular 
lines (see middle panels of Eig. [1]). 

The SSI variability on the timescale of solar rotation is 
mainly caused by transits of magnetic features across the visible 
solar disc as the Sun rotates. The spot and facular contributions 
to the irradiance rarely cancel each other at any specific moment 
of time. Hence, in contrast to the case of the 11-year variabil¬ 
ity, the amplitude of the irradiance variability on the rotational 
timescale does not drop in the visible part of the spectrum and 
is still signihcant in the red and IR spectral domains (where it is 
mainly associated with umbral contributions). Consequently the 
role of the wavelengths affected by the molecular lines is signif¬ 
icantly smaller than in the case of the 11-year variability. 

In our model the spectral prohle of the irradiance variabil¬ 
ity on the centennial timescale is directly given by the spectral 
prohle of the facular brightness contrast (because the same at¬ 
mospheric model is used for faculae and network, see Sect. 13- 
In the infrared the centennial SSI variability is larger than the 
11-year SSI variability (since there is no competing effect from 
spots) but is substantially smaller than the variability on the time 
scale of solar rotation. 


5. Effect of the Fraunhofer lines on spectral 
irradiance variability 

Despite the difference in overall shape, all three spectral prohles 
of the irradiance variability, introduced in Sect. 13 have similar 
spectral features which are mostly attributed to different molec¬ 
ular bands and strong atomic lines. To pinpoint the contributions 
of variations in molecular and atomic lines to these prohles we 
calculated intensity spectra (A, r) by hrst putting the opacity 
in molecular lines to zero and then the opacity in all Eraunhofer 
lines to zero. Eor each of the computations Eq. ([T]l was then used 
to produce “SATIRE-NESSY no molecules” and “SATIRE-NESSY 
no atoms” SSI time Series, respectively. These SSI time series 
were then substituted into Eqs. (|3-(|3 to recalculate the spec¬ 
tral prohles of the irradiance variability as they would be mea¬ 
sured without contributions from molecular/Eraunhofer lines. To 
maintain a consistent normalisation of the spectral prohles, the 
TSI time series used in Eqs. (|2ll-(|3 was calculated taking the 
molecular and atomic lines into account, i.e. is always the same. 

5.1. Molecular lines 

In Eig. [3] we present a comparison of the spectral prohles 
of the SSI variability calculated employing the “SATIRE- 
NESSY” (produced with the full linelist) and “SATIRE-NESSY 
no molecules” SSI time Series. In line with the discussion in 
Sect. I3.11 molecular lines mainly affect the near UV, violet, and 
blue spectral domains. They almost double the SSI variability 
on the 11-year and centennial timescales between 300 and 450 
nm. As mentioned in Sect. |3 this spectral domain also happens 
to be the major contributor to the TSI variability on the 11-year 
timescale. By integrating the spectral prohles of the centennial 
and 11-year variability over the wavelengths one can estimate 
the contribution of the molecular lines to the TSI variability. Our 
calculations indicate that 23% of the TSI variability on the 11- 
year timescale is attributed to molecular lines. On the centennial 
timescale (bottom panel of Eig. [3 the relative contribution of 
the 300-450 nm spectral domain to the TSI variability is smaller 
than in the case of the 11-year variability. Consequently the con- 
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TIMESCALE OE SOLAR ROTATION 



ACTIVITY CYCLE TIMESCALE 



CENTENNIAL TIMESCALE 



Wavelength [nm] 


Fig. 3. The spectral profiles of the irradiance variability calcu¬ 
lated on the solar rotational (upper panel), 11-year activity cycle 
(middle panel), and centennial (lower panel) timescales. Plotted 
are the values calculated with “SATIRE-NESSY” model taking 
the full linelist into account (red curves) as well as putting the 
opacity in molecular lines to zero (black curves). 


tribution of molecular lines to the TSI variability on the centen¬ 
nial timescale is also slightly lower, being equal to 15%. 

In contrast to the 11-year and centennial variabilities, molec¬ 
ular lines do not play an essential role in the variability on the 
timescale of solar rotation. This is mainly because of the spot 
contribution which keeps the rotational variability high even in 
the absence of molecular lines. 

Molecular lines are also responsible for the strongest peak 
in the SSI variability on all three timescales considered in this 
study, namely the 380-390 nm peak which is linked to the CN 
violet system. The second strongest peak in the SSI variability 
on the 11-year and centennial timescales is associated with the 
CH G-band. 


5.2. Fraunhofer lines 

In Fig.|4]we compare the spectral profiles calculated employing 
the “SATIRE-NESSY” and “SATIRE-NESSY „„ SSI time 
series. The elimination of the Fraunhofer lines has a profound 
effect on all three spectral profiles up to 600 nm. The ionisation 
edges clearly visible in the brightness spectra and contrasts in 
the right panels of Fig. [1] are also prominent features of spectral 
profiles of continuum variability. 

Notably our calculations indicate that the continuum SSI 
variability in the UV and visible is in antiphase with solar activ¬ 
ity o n the 11 -yea r and c entennial timescales. This is consistent 
with lUnruh et al] (Il999h . who also reported an antiphase con¬ 
tinuum variability on the 11-year timescale for the 300 -500 nm 
spectral domain (cf. also lCriscuoli & lJitenbroel3l2014ll . 

Furthermore, since the 11-year and long-term TSI variability 
calculated without accounting for spectral lines is out-of-phase 
with the solar cycle, the contributions of Fraunhofer lines to the 
TSI variability on these timescales is larger than 100%, being 
250% and 140% for the 11-year and long-term timescales, re¬ 
spectively. By the same token the increase of the TSI at maxi¬ 
mum of the activity cycle compared with its minimum is directly 
attributed to the variability in spectral lines. We note, however, 
that the exact behaviour of the continuum variability strongly 
depends on the temperature structures of the quiet Sun and mag¬ 
netic features in the deepest photospheric layers. These layers 
provide a very small contribution to the emergent line spectra, 
so that it is difficult to reliably constrain their temperature struc¬ 
tures in ID semi-empirical modelling, especially taking intro 
account the uncertainties in the measurements of the IR solar 
irradiance which emanates from the deep photospheric layer s 
(iThuillier et alJ 120141 : iBolsee et al.ll2014t iThuillier et al.ll2015h . 
Consequently our model does not allow us to unambiguously 
rule out an in-phase continuum variability. At the same time our 
calculations clearly indicate that the variations in spectral lines 
govern the TSI variability on all three considered timescales. 


6. Discussion and Conciusions 

We have estimated the contribution of the Fraunhofer lines to the 
solar irradiance variability driven by the evolution of the surface 
magnetic field on the solar rotational, 11-year activity cycle, and 
centennial timescales. Our calculations indicate that the solar ir¬ 
radiance variability in the UV, violet, blue, and green spectral 
domains is fully controlled by the Fraunhofer lines. The highest 
peak in absolute SSI variability on all considered timescales is 
associated with the CN violet system between 380 and 390 nm. 
Furthermore, the molecular lines strongly enhance the SSI vari¬ 
ability on the 11-year and centennial timescales. For example, 
according to our model, almost a quarter of the TSI variability 
on the 11-year timescale originates in molecular lines. 

Such a substantial contribution of molecular lines to the ir¬ 
radiance variability is due to the strong sensitivity of molecular 
concentrations to temperature changes and consequently an en¬ 
hanced molecular depletion in hot facular regions compared to 
the colder surrounding regions of the quiet Sun. This is a well 
known effect, which leads to the appearance of bright points in 
molecular lines at the locations of magnet ic flux concentrations 
(ISteiner et alj|200il ISchtissler et al.ll2()03h and makes molecu¬ 
lar lines a sensitive t ool for mapping photosp heric temperatures 
and magnetic fields (iBerdvugina et al.N2005h . However, to our 
knowledge, it is the first time that molecular lines are directly 
linked to the enhancement of the irradiance variability. 
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TIMESCALE OE SOLAR ROTATION 



ACTIVITY CYCLE TIMESCALE 



CENTENNIAL TIMESCALE 



Fig. 4. The same as Fig.[3]but now the black curves represent the 
values calculated putting the opacity in all Fraunhofer lines to 
zero. 


Some words of caution are needed here since our model 
is based on the NLTE calculations with static ID atmospheric 
stratifications of the quiet Sun and magnetic features. Such 
calculations are still being actively devel oped and are getting 
highly sophisticated (see e.g. discussion in iRutten & Uitenbroek 
2012). At the same time they are dee med (e.g. iKoesterke et ~ 


2008 ; lUitenbroek & Criscuolil 1201 ih to be inherently unreli¬ 


able for diagnostics o f solar atmospheric properties (see how¬ 


ever 


Vitas et 


(IStenflo et al 


LI|2009Ii . unless differential techniques are used 


19981) . The ID models cannot directly account 


for the hot walls of flux tubes which to large ex tent determin e 
the centre-to-limb variation in facular brightness (ISteineill200-5h . 
Furthermore, molecular concentrati ons depend on the temp er- 
ature in a non-linear way (see e.g. [B^vugina et alj|2003l) so 
that their response to horizontal temperature fluctuations can¬ 
not be properly averaged out in ID calculations. Also, for some 
molecule s the assumption of instantan eous chemical equilibrium 
may fail dAsensio Ramos et al.ll2003h . 


The temperature structures of ID solar atmospheric mod¬ 
els are specifically adjusted to satisfy as many observational 
constrains as possible. This makes ID models quite success¬ 
ful at reproducing the obser vations relevant for calculations of 
the SSI variability (see e.g. lUnruh et"^ Il999t iFontenla et al.l 
1201 ll and references therein). Additionally, despite their huge 
potential, self-consistent 3D radiative transfer calculations of 
the entire solar spectrum based on the results of MHD simu¬ 
lations and capable of reproducing h igh-resolution obse rvations 
are not yet readily available (see e.g. lAfram et alJl201 ih . Hence, 
despite the obvious limitations, the semi-empirical ID calcu¬ 
lations still remain at the forefront of irradiance modelling. At 
the same time 3D MHD simulations have been gradually reach¬ 
ing a new level of realism and there are also inter mediate ap¬ 
proaches in between ID and 3D radiativ e transfer dAvres et al.l 
l2006t iHolzreuter & Solankill2012l l2013h . It would be interest¬ 
ing to repeat the exercise presented in this paper when such cal¬ 
culations are applied to modelling of the solar irradiance vari- 
ability. A study along t hese l ines was recently performed by 
ICriscuoli & UitenbroekI d2014ll . who, however, only calculated 
the facular brightness at continuum wavelengths. 

Our result indicates that proper calculations of opacity in the 
solar atmosphere are of high importance for modelling of the 
solar irradiance variability. The inclusion of Fraunhofer lines is 
absolutely crucial for reproducing variations in TSI and SSI at all 
timescales from the solar rotation period to centuries. For quan¬ 
titatively accurate computations molecular lines also need to be 
taken into account, in particular CN and CH. 
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